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ABSTRACT 

Many workers have found that the recollapse of a dark matter halo after 
decoupling has a self-similar dynamical phase. This behaviour is maintained 
strictly so long as the infall continues but it appears to evolve smoothly into 
the virialized steady state and to transmit some of its properties intact. The 
density profiles established in this phase are all close to the isothermal inverse 
square law however, which is steeper than the predictions of some n-body 
simulations for the central regions of the halo, which are in turn steeper than 
the density profiles observed in the central regions of some galaxies; particu- 
larly dwarfs and low surface brightness galaxies. The outer regions of galaxies 
both as observed and as simulated have density profiles steeper than the self- 
similar profile. Nevertheless there appears to be an intermediate region in 
most galaxies in which the inverse square behaviour is a good description. 
The outer deviations can be explained plausibly in terms of the transition 
from a self-gravitating extended halo to a Keplerian flow onto a dominant 
central mass (the isothermal distribution can not be complete) , but the inner 
deviations are more problematic. Rather than attack this question directly, 
we use in this paper a novel coarse-graining technique combined with a shell 
code to establish both the distribution function associated with the self-similar 
density profile and the nature of the possible deviations in the central regions. 
In spherical symmetry we find that both in the case of purely radial orbits 
and in the case of orbits with non-zero angular momentum the self-similar 
density profile should flatten progressively near the centre of the system. The 
NFW limit of —1 seems possible. In a section aimed at demonstrating our 
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technique for a spherically symmetric steady state, we argue that a Gaussian 
distribution function is the best approximation near the centre of the system. 

Key words: dark matter: galaxy formation : distribution function: coarse 
graining 



1 INTRODUCTION 

The study of the radial infall of dark matter and the realization that the evolution becomes 
self-similar has by now a long history (see for example QHenriksen and Widrow 1999|) for 



a summary of the relevant references and a nearly state of the art statement of our under- 
standing ). The salient theoretical behaviour is well known if not completely understood, 
but it yields density profiles that are theoretically distinct from the inverse square law but 
observationally indistinguishable from this law. However the NFW law from n-body simu- 
lations ([Navarro, Frenk and White 1996|) even if slightly steeper when resolution effects are 



taken into account (e.g. ( Moore 98 ); but see also ( [Kravtsov 98|) ) is distinctly flatter than the 



inverse square law in the central regions. Moreover in some galaxies (low surface brightness 
or LSB galaxies, DE Blok et al. 2001; some dwarf ellipticals, Stil 1999; and even brighter 
elliptical galaxies compared to the fainter objects, Merritt and Cruz 2001) the observed 
density profile is significantly flatter in these regions than even the n-body profiles (except 
( |Kravtsov 98j )). Both observed and simulated profiles are steeper than the inverse square 
law in the outer regions, but this can be understood ( |Henriksen and Widrow 1999| ) in terms 
of secondary accretion after most of the halo mass has fallen in. 

We take the view in this paper that if finite resolution effects are at work in the centre 
of the simulations, and perhaps also in reality since obviously real particles can not yield an 
infinite density cusp; then we may gain some insight by solving analytically the Collisionless 
Boltzmann Equation (CBE) for the distribution function (DF) in an series expansion in 
powers of the inverse 'smoothing length' (to be defined precisely below). In such an expansion 
the lowest order is the coarsest grained approximation and higher orders yield progressively 
finer grained information. It transpires that in this way we can deduce the radial distribution 
function necessary to maintain the strict self-similarity and we may place some constraints 
on modifications to this DF that can produce flattening. 

Our method of coarse graining is by non-canonical transformation on the phase space 
(we make a special choice of a stretching transformation here, but the method is more 
general) which produces a definite equation to be satisfied by each term in our series. We il- 
lustrate the method first by applying a stretching transformation to a spherically symmetric 
isotropic equilibrium. The bulk of the paper is devoted to exploring the self-similar 'equi- 
librium'maintained by the continuing accretion onto a relaxed core ( |Fillmore and Goldreich 
1984Q ; ( [Bertschinger 85| ); ( |Henriksen and Widrow 1999| ) 



A consideration of the entropy for the system allows us to suggest that the central DF 
in spherical equilibrium is Gaussian in velocity with a separable radial dependence, while in 
exact self-similar secondary infall the DF is everywhere proportional to the square root of 
the particle binding energy. Moreover the requirement that the coarse grained entropy be 
maximal predicts perturbations that flatten the inner regions of steep self-similar profiles 
and steepen the outer regions of flat self-similar profiles. We show also that the addition of 
angular momentum to spherically symmetric orbits does not change these conclusions. 

* henriksn@astro.queensu.ca 
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2 COARSE GRAINING OF THE CBE BY NON-CANONICAL PHASE 
SPACE TRANSFORMATION 

We wish to coarse-grain the CBE together with the mean field Poisson equation in some 
systematic fashion. A simple smoothing of phase space to get a coarser grid over which to 
define the DF suffers from non-commutation with the operators in the CBE and leads to 
a coarse grained DF that is defined only qualitatively (e.g. QBinney and Tremaine 19871) ). 



Our idea is to use a coordinate transformation on phase space that does not preserve phase 
space volume and so is non-canonical. Although this can be applied to evolving systems it is 
not very well defined because of a lack of information concerning the initial conditions of the 
coarse grained (time averaged for an evolving system) DF. Thus we will restrict ourselves 
to systems that have attained 'equilibrium' in some appropriate variables (see below). 

The simplest example of this idea is to apply a stretching transformation to phase space, 
and in fact this is the only transformation that we have explored (the scale factor depends 
on time in the self-similar examples). However more general non-canonical transformations 
may be found that are useful so we continue with this label. 

This can be readily carried out on an equilibrium spherical collisionless system for which 
the governing equations are; 

Vd r f - d r <S>d v f = 0, (1) 

and 

«9 r (r 2 «9 r $) = AnGp r 2 = AnGr 2 J f Anv 2 dv, (2) 

where / is the normal phase space mass density and v is the magnitude of the isotropic 
velocity. A scaling or stretching transformation can be applied to the phase space as 

r <— rt v <— uv, (3) 
f<-ffi $ <- 4nGp£ 2 u 3 $, (4) 
p <- pu 3 p, (5) 

where £, p and u represent constant factors that we may take here to be dimensionless 
numbers. The phase space volume thereby transforms by the factor £ 3 u 3 and we will take p 
equal to the reciprocal of this volume times a fiducial mass m in order to conserve mass (for 
unbounded systems). Thus £ and u are the theoretical smoothing 'lengths' on phase space. 
Finally we scale / according to / <— // (4nG). 

The equations are the same in the stretched variables except that the CBE becomes 

vd r f ~ K^drQdJ = 0, (6) 

where 

K = u 2 £/4nGm. (7) 

As either or both of the phase space variables are stretched this parameter becomes large. 
We coarse grain the system by assuming a convergent expansion of the form 

f = T, i= ,f t {r )V )Tl-\ (8) 

so that the lowest order contains the least information about the DF. A similar expansion 
follows for the density and the potential. 

Proceeding in this fashion we find first that f = f (v 2 ) and consequently that 

po = J fo^v 2 dv. (9) 

The limits in this equation would normally be set so as to include just the bound particles. 
However the stretched specific energy (E <— u 2 E) is related to the other stretched quantities 
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as E = v 2 /2 + 7?. _1 $ and so is purely kinetic to lowest order. Thus we take the upper limit 
for p to be simply v m , some maximum speed at r = 0, and the lower limit to be zero. 
The integral so defined is denoted p {v m ). Consequently one finds the harmonic potential in 
lowest order (so that all particles are bound and turn where E = 1Z~ l Q) 

$o = p r 2 /6. (10) 
The first order term in the DF expansion is now found to be 

and consequently 

4?r „ ji r Vm ..2 # 

dv 

*i = -p r J o v—dv, (13) 

assuming no central point mass. The next order terms in these series may be found to be 
proportional to r 4 , since for example 

2 P yd 2 f 2 df 

drh = —dW + 5 pir ^ (14) 
and so on. In this way we are generating an expansion of the type necessary in general for 
polytropes (e.g. ( |Chandrasekhar 1957|) p. 94), but here we are able to relate the coefficients 



Pi = —Par /_ v —dv, (12) 



directly to the DF that is not necessarily a polytrope. In essence since most of the physical 
system will be at small r for an arbitrarily stretched coordinate, we are able to re-interpret 
the polytropic expansions in terms of progressively finer graining. 

The second order term in the DF is proportional to r 4 , but more importantly it involves 
the second derivative of f with respect to v 2 . This behaviour will continue to all orders where 
the highest order derivative of f appearing is equal to the coarse-grained order. One must 
therefore be careful to ensure that these derivatives are sufficiently well-behaved to all orders 
that a divergent density is avoided, which would render the expansion invalid. For example if 
we try to imitate a polytrope by setting (A is a positive constant) f„ = Aty+v^ — v 2 p n ~ 3 / 2 \ 
then all derivatives are well-behaved for v < v m at the cost of a discontinuity in the DF at 
v m . The quantity ip appears as an negative constant potential (which is stretched by the 
factor u 2 ) to be added to $ D and in principle we can keep only the negative energy particles 



by setting v m (r) = \j2(vp — 1Z 1 p Q r 2 /6). However the presence of 1Z in this expression breaks 
the ordering of our expansion and although it may improve the rapidity of convergence, we 
continue to treat constant and leave the corrections to higher orders. 

To obtain a continuous transition to zero at v m , we might take for f a one of: a Gaussian 

DF 

f a = A>™-^/ 2<7 ; (15) 
or a King type DF (a constant potential is absorbed into v 2 ^) 

f = A(e^- v2)/2a - 1). (16) 

Ultimately the choice of f Q must be that which best reflects the coarse-grained image 
of the system in question. This seems to require us to know something about the gross 
properties of the system 'a priori' so that in fact our procedure is more like a progressive 
'fine graining'. Each choice of coarse-grained DF f a will generate by this procedure (provided 
that it may be carried out properly as discussed above), a progressively more precise DF as 
the fine graining terms are added. 

We may seek general principles concerning the form of f a however. We note that the 
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choice of the Gaussian is unique in the present context in that it will maintain this de- 
pendence in the DF to all orders. In addition it will maximize a version of the statistical 
entropy flNakamura 200QO ,( [Lynden-Bell 1967] ) . The expansion is well defined to all orders 



for the Gaussian, but we can only expect it to apply to the central regions of finite systems 
(otherwise one will not conserve mass as the scale is stretched). 

The expansion of the Gibbs' statistical entropy (essentially the generalized 'H function' 
of Boltzmann when there are no correlations) in the form (f is normalized by the total mass 
of the system) 

S = J dr/ lnl// -^(l-lnl// ) + -.., (17) 

shows that since f Q is a small number (non-degenerate system) and since fi < ( actually 
true for all three cases mentioned above although we require n > 3/2 for the polytrope), 
the entropy decreases as the fine graining terms are added. This seems as it should be. 
Moreover with the choice of the Gaussian for f a the DF is separable. This makes the velocity 
distribution the same at all r as one might expect in the central equilibrium region. 

Our conclusion in this section is that the best behaved coarse graining expansion is 
consistent with the proposition that a Gaussian should be the DF in the central regions 
of finite spherical systems that have undergone relaxation and hence the evolution towards 
maximum entropy. Our arguments do not of course prove this proposition, but they support 
the statistical arguments ([Nakamura 2000|) , ([ShuT978|) , (|Shu 1987|) , ( [Lynden-Bell 19671 ). We 



have essentially used the idea that in a relaxed system we expect the fine grained DF and 
the coarse grained DF to be the same. That is that the DF should be independent of cell 
size. In fact the principal advantage of our method is that the coarse graining expansion 
dispenses with the necessity in the combinatorial statistical treatment of making an arbitrary 



distinction between microcells and macrocells ( Shu 1978 ). Moreover it is in agreement with 
the maximum entropy approach of Nakamura ( Nakamura 2000 ). 

It may seem contradictory to infer a DF that remains the same at all radii as found 
above while also restricting our proposition to the central regions of the system. In fact in 
a completely relaxed system the DF should be the same at all radii as in fact is true for 
the Gaussian DF. But we know this to yield an infinite system ([Lyndcn-BcII 1967]) , and 
so presumably deviations from the relaxed state increase with radius in any finite system. 
This means that at large radii the coarse grained DF and the fine grained DF no longer 
coincide and the coarse Grained DF need no longer satisfy Liouville's theorem. Moreover 
our stretching transformation can not be applied too close to the boundary of a finite system. 

These ideas suggest for example that polytropes are not in a maximum entropy condition 
but rather have to be constructed rather carefully. The density profile consistent with a 
Gaussian DF does not have a central cusp so that simulations that predict such cusps 
are presumably not sufficiently relaxed near the centre. Our investigations of the next two 
sections suggest in fact that progressively finer resolution of the central regions of self-similar 
infall produce flatter density profiles. 



3 COARSE GRAINING OF SELF-SIMILAR INFALL WITH RADIAL 
ORBITS 

When we restrict ourselves mainly to radial orbits, we will introduce as usual the canonical 
distribution function F(t, r, v r ) such that the complete phase space mass density is given by 

/ (4 7 r 2 G)- 1 <5(i 2 ) F, (18) 

where j 2 = r 2 (vg + f^). The mean field treatment is then given by the combined CBE and 
Poisson equation in the forms(see e.g. Henriksen and Widrow, ibid) 
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d t F + v r d r F -d r $d Vr F = 0, (19) 
and 

d r (r 2 d r $) = J F dv r = 4ttGV 2 p, (20) 

where $(t, r) is the gravitational potential and p(t,r) is the mass density. 

In the previous section we pointed out that our procedure made the most sense when 
applied to a system that is in 'equilibrium' since otherwise the coarse graining would also be 
over the history of the system. However a time-dependent case which is nevertheless in a kind 
of equilibrium is the self-similar phase of the formation of the system. In fact by changing to 
variables appropriately scaled b y the time (|Fillmorc and Goldreich 1984|) ; (|Bertschinger 85|) ; 



( Henriksen and Widrow 1997| );( |H"enriksen and Widrow 1999 ) the equation is transformed to 
an equivalent steady problem. 

Here we adopt these variables in a formalism given by Carter and Henriksen (1991; see 
also QHenriksen and Widrow 1995|) ; ( Henriksen and Widrow 1999|) ; (Henriksen 1997|) ) where 



we adopt an exponential time T 

e aT = at, dT/dt = e~ aT . (21) 

It is convenient to take time to be measured in units of some fiducial value so that a, the 
stretching scale for time, is also a dimensionless number. We introduce the scaled phase 
space variables X, and Y for the radius and radial velocity respectively as (in the previous 
section the scaled quantities were not given new labels but we continue to do this here to 
agree with previous practice) 

X = re ~ (<5/a)aT , Y = v r e-( s/a ~ 1)aT , (22) 

where 5 provides for an independent spatial stretching scale. However in a self-similar system 
this ratio is constant. In addition we scale the radial canonical phase space DF according to 

P(X,Y) = Fe- (s/a - 1)aT , (23) 

where the lack of a dependence on T in the scaled DF P is the condition for self-similarity. 
The potential is correspondingly scaled to ^(X) where 

* ee e ~ 2(<5/a ~ 1)aT $, (24) 
and the density is scaled to 0(X) where 

9 ee pe 2aT . (25) 
Consequently 

6 = — '— / P dY. (26) 



AnX 2 „ 

We observe that the phase space volume element ArAv r has been scaled according to 
ArAv r = AXAYe^ a ^ aT . (27) 
The scaling 'preserves' mass since 

PAXAYe {3S/ - 2a)T ee FArAv r = Am, (28) 

together with 35 — 2a = fi where \x is the mass scale (e.g. ( [Henriksen and Widrow 1995) )) 
give the correct time dependence for the mass element Am during the self-similar phase 
(which is only constant for Keplerian self-similarity wherein 5/a = 2/3). 
Equations (H) and © now become respectively 

(*-l)P + (^-(VW - ((*-l)r + i§W = 0, (29) 
a \ a a J \ a a dX I 
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P alY. (30) 

These are the equations to be coarse grained in this case. 

We make use of the parameter a in the transformation to the self-similar phase space 
variables to effect the coarse graining. We note that equation (p7|) shows that a fixed volume 
element AX AY at a fixed time corresponds to an ever larger volume element Af r Ar as a 
is increased while holding the ratio 5 jot constant so as to maintain the similarity class. This 
is true so long as the similarity class is > 1/2, which is usually the case of interest since the 
Keplerian value of 2/3 tends to be the minimum value encountered. The value 5 /a = 1/2 




gives a canonical transformation by equation (|27j) and no change in the volume of phase 
space is effected. Increasing a amounts to a parametric way of stretching time and space 
so that it becomes the theoretical smoothing length parameter. Normally we take a = 1 
which can be regarded as the fine grained limit. Here however we use the freedom in its 
absolute value to allow it to adopt large values and thus yield a coarse grained limit (with 
the similarity class fixed the spatial scale is also being stretched in proportion). 
The expansion we use is again of the form 

P(X,Y) = i: i= oP l a-\ (31) 

and similar expansions apply to the density and the potential. We begin with the zeroth 
order equations which become (it is easy to write the nth order formally but this is not very 
transparent) 

{--l)P -{-)Xd x P - (- - l)Yd Y P = 0, (32) 
a a a 

P D alY. (33) 

We readily find P Q by the method of characteristics to be 
P = P 00 (C)X {1 ~ a/s \ (34) 

where P oa is an arbitrary function of (, which is constant on characteristic curves of P a . 
These curves are in turn given by 

Y 

C = X (l-a/S) ' ( 35 ) 

and the actual arc-length s of a characteristic in phase space may be taken to be 

s = ^lnX, (36) 
o 

and this can be used to give X(s), Y(s) if desired. 

If we recall equation (|26|) and combine it with ([35]) for dY and (|34"D for P then we may 
write 

0o = J Poo(C) d(. (37) 

Thus the zeroth order coarse graining already produces the self-similar density profile in r, 
namely r~ 2a l 5 (cf ( Henriksen and Widrow 1995|) ; ( Henriksen and Widrow 1999|) ). This is not 
surprising once the self-similarity is imposed since it is essentially a dimensional argument 
(e.g. ([Henriksen and Widrow 1995 )). Since however at this level no potential enters the 



calculation, it does suggest that such elements of a simulation as smoothing length may not 
be essential to finding this profile if the other elements that impose the self-similarity (such 
as initial density profile) are in place. We shall see in a shell code application below that 
there is some evidence for this. 
The zeroth order potential is 
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J v(2-2a/S) 

^-(3-f)(2-f)- lX ' (38) 



where I OQ is the integral over P oa occuring in equation (|37|) and 7 is defined to be positive 
when a 1 5 > 1. The cases a/5 = 1 and a/<5 = 3/2 are logarithmic rather than power law and 
must be treated separately. The potential is not in general simply harmonic (except in the 
steady case where a = that does not concern us here. An expansion in positive powers of 
a would allow a 'fine graining' about this state.) 

Finally at this order we note that the energy of a particle scales as E = e 2 ^~ a ^ T e where 

e = Y 2 /2 + ty(X), (39) 
and with e D = e 00 X 2 ( 1 ~ a / <5 ) we find that £ is related to the zeroth order energy through 



C = ± v /2(e 00 + 7 ). (40) 



Consequently the limits of the integration over £ at this order are [— +y/%y\ when 7 
is positive in order to include only the bound particles, while when 7 < we have a case 
similar to the harmonic potential of the previous section and ( varies between plus or minus 
some maximum value. 

We now proceed to the first order terms for which the following equations must be solved 

(6/a-l)P 1 --Xd x P 1 -(6/a-l)Yd Y P 1 = -Y0 x P o + ^d Y P , (41) 
a aJi 



(42) 



d 

dX 

Once again the method of characteristics yields a solution as 

Pi = x (i-2a/ S) {{a/5 _ 1)c(Poo _ cp ,j + Ioo pJ {3 _ 2a/ 5)) = PniOX^ 1 - 2 ^, (43) 

where the shape of the characteristics (for Pi) is the same as for P a above and the prime 
indicates differentiation with respect to (. In addition the form for 8\ follows as 

v(-3a/s) r 

0i = J Pn(0 d(- (44) 

Thus the solution to first order in the coarse graining parameter (i.e. a finer grained solution) 
takes the form 

P = X^- a ^(P oo (0 + a- 1 P 11 (()X-^), (45) 

v-(2a/8) 

9 = -——^I 00 + a - 1 X- a / 5 I 11 ), (46) 

where In is the integral occurring in 9\, which is not however positive definite in this case. 

Now what is remarkable about these forms is that they suggest that at small X (large 
t and/or small r) there are deviations from the coarse grained self-similar behaviour as 
we increase the resolution or information content. Indeed should P\\ be negative then there 
would be a flattening of the density profile in this limit. However we know from the numerical 
simulations that there is a region where the self-similar behaviour is maintained so long as 
the infall continues. This leads us to ask a question similar to that which was posed in the 
previous section. There we asked for separability in the DF so as to yield the same velocity 
distribution at all radii. Here we can ask for the DF which maintains the self-similar density 
profile for all X. This requires Pn = since then all higher order terns in the expansion will 
vanish and the density profile is uniquely self-similar. We can write this term from equations 
(M) and (HUD in the form 
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Pu = ia/ S- lK P^-2 d J^), (47) 

so that the condition that all higher order terms vanish is simply 

P QO = const x \e 00 \ 1/2 = const x \( 2 /2 - 7| 1/2 . (48) 

Retracing the various definitions and scalings, we can discover that this last result gives 
for the canonical radial DF 

F(E) = const x \E\ 1/2 , (49) 
and hence the full DF is 

/ = const, x (4tt 2 ) _1 5(f)\E\ 1/2 . (50) 

This tends to confirm the conjecture of ( [Henriksen and Widrow 1999| ). This was based 
largely on the argument for a continu ous transition to the steady state (a steady state found 
also in ( [Henriksen and Widrow 1995|) ) together with some weak evidence from a shell code 
simulation. Our present argument is related since in the limit that a — > oo we obtain a 
steady state, and in fact it is readily verified from (|35|) that ( is independent of time. By 
requiring all higher orders to vanish we have effectively selected this steady state. 

In the next section we shall show some increased numerical evidence for this DF. The 
evidence is best when a point mass is allowed to be present at the centre, even if it is a 
negligible fraction of the total halo mass. This might be termed a 'black hole' but in reality 
it simply imitates the density singularity expected from the self-similar density profile. That 
profile does not yield a finite mass at the centre but it does predict a cusp that the numerical 
work has difficulty expressing in the absence of a central point mass. When the point mass is 
small compared to the halo mass it serves to imitate a cusp without a finite mass singularity. 
Thus we might well expect the DF in this case to be closest to the (— E) 1 ^ 2 behaviour (see 
below) . 

The real significance of this result is to show that deviations from the self-similar infall 
profile of radial orbits in spherical symmetry (a subsequent section will deal with more 
general orbits) will be expressed as deviations from the (— E) 1 / 2 law. These may well be due 
to incompleteness of the self-similarity at the centre, to angular momentum or to geometric 
effects. We have already seen that an abrupt cut-off at high binding energies does create a 
deviation in the density profile, but only to logarithmic order ( |Henriksen and Widrow 1999| ) 



More global changes must occur in the DF in order to produce substantial changes in the 
density profile even near the centre. 

Finally it is also of interest to consider the entropy expansion. We follow custom (e.g. 
( |Binney and Tremaine 19"87| )) in writing the Gibbs' entropy as 

S = -J flnfdr, (51) 

where dr = Airr 2 dr 2irjdj dv r /r 2 . This is the generalized (i.e. to inhomogeneous systems) H 
function of Boltzmann in the absence of particle-particle correlations. Then if we scale the 
entropy and mass of the system according to 

S = Se-^- 2a ^ T } (52) 

M = Me-( 35 " 2Q ) T , (53) 

we may find that 



S = M(\n(A7r 2 ) + (l-5/a) aT)- J dX dY PlnP, (54) 
where the scaled mass is related to the scaled density in the natural way 
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M = J AttX 2 6 dX. (55) 

Thus the only time dependence in this self-similar entropy is the explicit logarithmic increase 
with t (through the term in T) that is not very significant (it is in fact strictly zero when 
5/a = 1). The true entropy is increasing in direct proportion to the mass as is seen from 
the respective scalings. 

Let us now substitute the coarse graining expansion from equation (f46|) to find to first 
order 



S = M(\nAn 2 + (l-5/a)aT) + 

J dX X 2(1 " a/5) d( P 00 (In \jP QO + {a/ 5 - 1) lnX) 

+ / dX X 2{1 ~ a ^ d( ^lx~ a/s ((a/ 5 - 1) In X + In l/P OQ - 1) . (56) 
J a 

We consider first the undisturbed self-similar state wherein Pu = 0. The integrals over 
( and X are independent (when a/5 > 1 the integral over ( goes from —y/2r/ to +y/2j, 
when a j 5 < 1 it varies from an arbitrary negative number to the same positive number) 
so we may treat them as products. We observe that when a/5 > 1 the mixing part of the 
entropy (i.e. the integral over phase space) has a well defined value provided that a/ 5 < 3/2. 
However this mixing entropy should be positive and since the integral over X is dominated 
by the behaviour at small X when a/5 > 1 we see that this requires roughly that 

X^ s -V > exp - (| d( P 00 In \/P 00 )/ J d( P 00 ) . (57) 

In other words, as is usual we expect that there will be an inner limit to the extent of the 
steep self-similar state after which distortions must arise. 

When a/5 < 1 the negative mixing entropy will arise at large X so that the same 
inequality (|57|) yields an upper limit to the extent of the self-similar state. In general this 
limit may overlap with the non-self-similar starting conditions for the accretion of a dark 
matter halo. 

In the limiting case with a/5 = 1 (which by equations (|26"D and (^) has the density 
profile of the singular isothermal sphere), we see that the self-similar entropy is strictly 
constant in the scaled variables. This suggests a certain stability for this profile, but if we 
look at the coarse-grained mixing entropy we see that it diverges for an infinite system. 
Consequently such a state can never exist over all space for this reason as well of course 
because of the infinite mass that would imply. It does suggest that this might be the most 
stable self-similar state in a finite relaxed region however ( [Lynden-Bell 1967| ). 

We turn now to the fine graining term involving P n . This term represents a departure 
from the strict self-similar density profile as we have seen. Given however that it represents a 
finer graining of the system, we should expect its contribution to the entropy to be negative. 
Assuming that we are at small enough X (the fiducial X may be taken as the boundary of 
the core and the integral over X in this case is dominated by values at small X) that the 
bracket in the first order term of equation (^) is negative, then we require In > for a 
negative first order term. By equation ( pETD if P QO is symmetric in ( (such as is a power law in 
e oc ) then a sufficient condition is that the perturbed P OQ satisfy P ao oc |e 00 |^ n ^ where n < 1/2. 
This last condition together with equations (|37|) and (|38D shows however that the disturbed 
coarse-grained density will be flattened (0 O oc X~ 2a l s X( n ~ l l 2 ^ 1 ~ a l & * ) ) when n < 1/2. Thus 
the deviation from strict self- similarity will be such as to flatten the steep self-similar density 
profile from the centre outwards. One expects the disturbances to arise near the centre of 
the system by our considerations of the coarse-grained entropy. 
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For a/5 < 1, the integral over X is dominated by the behaviour at large X. Thus the 
requirement for a negative first order entropy (assuming the term in lnX to be dominant 
and negative) is now In < 0, and hence following the argument above and referring to 
equation (f|7|) we see that once again n > 1/2. However an examination of the dependence 
of o on X as above now shows that the profile will be steepened. In this case we can 
expect disturbances to arise at large X by our coarse-grained entropy argument. That is the 
deviation from strict self- similarity for the shallow initial profile is such as to steepen the 
profile from the outside inwards. This steepening is observed in the simulations to continue 
until the r~ 2 profile is obtained. 

Thus this exploration of the entropy function suggests that the flat self-similar behaviour 
(as determined by initial conditions) is unstable to the development of steeper behaviour 
, while the steep self-similar behaviour is unstable to flattening in the central regions. Ul- 
timately incompleteness dominates both at the centre and at the edge of the system as 
discussed above. Although this predicts no real 'universal' profile in the sense of (|Syer ancl 
White 1998|) , it does yield a universal qualitative behaviour in three sections: wherein the 
middle section is close to a profile of r~ 2 , an outer incomplete section is more like r~ 3 , and 
an inner section that is flatter than r~ 2 for a variety of reasons. This may be in accord with 
the simulations if not the observations. 

We should also observe that Evans and Collett ([Evans and Collett 1997|) have found a 
remarkable self-similar steady state solution to the collisional Boltzmann system that is an 
attractor and yields a flat r~ 4//3 cusp. This is suggested to be applicable to galaxy formation 
by way of mergers where the collisional ensemble is the set of merging clumps. Our approach 
focusses on the collisionless dark matter particles and stars and is therefore independent of 
this result. However it would be interesting to apply the present method of coarse graining 
to a collisional system to see how the Evans and Collett solution emerges. 

The analysis of this section has served primarily to test our coarse graining expansion 
against relatively well known results, although the results regarding the inevitable deviation 
from strict self-similarity are new. In the next section we turn to the more challenging 
problem of coarse graining the DF for a spherical system in self-similar evolution with 
velocity space anisotropy. That is j 2 ^ for each particle, but there is no net rotation of 
the system ( Henriksen and Widrow 1995|) . 



4 COARSE GRAINING OF SPHERICAL SELF-SIMILAR INFALL WITH 
ELLIPTICAL ORBITS 

In this section we show how our procedure may be applied to more complex systems. In 
principle we can dispense with spherical symmetry and consider axially symmetric and 
three dimensional systems. However the Green function solution of Poisson's equation then 
leads to algebraic complexities that tend to obscure the method somewhat. This will be 
attempted elsewhere. The present example is already new and physically interesting. 
The fundamental equations can be written in the form after Fujiwara ( |Fujiwara 1983| ) 

d t f + v r d r f+(£-d r <^d Vr f = 0, (58) 

d r (r 2 d r §) = Att 2 G J dj 2 J fdv r , (59) 
where the particle density is given by 

P = ^j d fj ( 60 ) 

We now define the anisotropic analogue of the self-similar radial infall of the previous 
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section. The definitions of X, Y and T are as given previously, while we introduce the extra 
scaled phase space variable Z according to 

Z = fe~ XT , (61) 
where on dimensional grounds we require 

A = 45 - 2a. (62) 

In addition we use the scaled potential \I/ as in equation (|24"D while the scaled DF is written 

as 

P(X, Y, Z) = An 2 Gfe^ s/a - 1)aT . (63) 
Then the equations that define an anisotropic self-similar infall model (ASSIM)are 

-(S--l)P+(---X)d x P- ((- - 1)Y + - A)) dyP-(d-2)Zd z P = 0, (64) 

a a a y a a dX X 6 j a 

and 

d ( #%\ = I iZ d YP. (65) 



dX \ dX / 

Moreover for the same scaling as in ( p5|) , the scaled density becomes 

9 = ^ J dZ dYP. (66) 

Once again the coarse graining consists in expanding all quantities as in equation fl3~T]). 
In this case we have scaled the phase space volume to 

ArAv r Af = AXAYAZe (G5/a ~ 3)aT (67) 

so that the condition for coarse graining at fixed T in the limit a —>■ oo remains 5/a > 1/2. 
One should note that should 5/a < 1/2, the expansion is about the fine grained limit since 
infinite a yields P Q as the exact result. 

Proceeding with the expansion the zeroth order equations become; 

AAA A 

(3- - l)P + -Xd x P + (- - l)Yd Y P + (4- - 2)Zd z P D = (68) 
a a a a 

U x ^)- i iziYp °- m 

We can solve for P Q by the method of characteristics to find 

Po = Poo(Ci,C 2 )e- (35/a - 1)s , (70) 
where the characteristic constants are 

Ci = Y/X {1 - a,&) (71) 

Cl = Z/X^- 2a ' 5 \ (72) 

and the arc-length along the characteristic s may be taken as in (|36|). 

It is interesting already to note that the scaled density may be written as 

6 = 7iX- 2a/5 I 00 (73) 

= nX- 2 */ 5 J P 00 (Ci,C 2 RidCf, (74) 

which shows that (provided the integral I oa goes over the same set of characteristics at every 
X) the usual self-similar density profile holds even to zeroth order, just as for the radial 
orbits. Since there is no potential term in the governing equations at this order, this must be 
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set solely by the initial conditions, which determine the similarity class a/5. Consequently, 
the zeroth order potential in this case is as in equation ( |3lf ) except that the integral I OQ is 
here defined as the integral over d and (2 that appears above in the scaled density. 
Using the same notation as in the radial case, the scaled energy is 

£ = * + T + 2^' (75) 

and so in zeroth order with e Q = e 00 X < - 2 ~ 2a//<5 ) as for radial orbits 

C 2 C 2 

e 00 = -l+ ( f + f, (76) 

where 7 is defined as before in terms of the current I ao . 

Proceeding to the first order in the expansion we find the governing equations to be 

^Xd x P 1 + (--l)Yd Y P 1 + 2(--l)Zd z P 1 = -(— -l)P a 
a a a a 




J dZ dYP v (78) 



Thus the form of the characteristics remain the same as in equations (71), (72) and the 
solution for Pi by the method of characteristics yields (plus a term that may be absorbed 
into the zeroth order since it has the same dependence on s) 

Pi = -e- 3|s Pn(Ci,C 2 ) 

= -X- 3 P u (Ci,C2), (79) 
where 

Pu = [ij ~ l)Ci 2 - (3 _ Io 2 ° a/5) + C 2 2 ) d (l P 00 + 2{- 5 - 2)Ulde 2 Poo - (3 - j)P oa . (80) 

Consequently the solution to first order in the coarse graining parameter for the DF is 

P = PooiCu&X-^V - -X- 3 P n (d,Ca), (81) 

a 

and the corresponding density profile becomes formally that of equation ( fl6|) but for a 
numerical factor, namely 

6 = nX- 2a ' s (l 00 - -I n X-^ s ) , (82) 



a 

Here however the integrals are defined as la = J Pad(id( 2 . 

We arrive then at a conclusion similar to the one found for purely radial orbits, namely 
that deviations from the self-similar density profile must arise as X — > 0. With Pu > 
for example so that In > 0, we see that the profile will be flattened near the centre of 
the system as the resolution increases. Should Pu = then we have the condition on P OQ 
that the system remain self-similar at all X to all orders. This requires redefining P QO to be 
Poo + ct X P \ + • • • where the P , are the contributions to the zeroth order variation from the 
jth orc [ er _ jf these functions are taken to be zero, the argument is exact. 

The argument based on entropy (that we do not reproduce here for reasons of brevity) 
suggests that there should be flattening at the centre and so we turn to consider the condition 
P\\ > 0. Just as in the radial case we may hope to establish the DF that maintains the 
coarse-grained self-similarity exactly and to put some constraints on deviations that maintain 
Pn>0. 
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The equation P n = is a partial differential equation for the exact self-similar DF. We 
may solve this equation once again by the method of characteristics to find 

,(3-a/W (83) 




-7j(C 2 2 )- 6 , (84) 
and setting v = £| for convenience 

± 2(a/S - 2)f = 1 (85) 

In these formulae the constants are 

a/6-1 
= 2(a/5-2)' 

7 = — . (86) 

The preceding solution holds for a/5 ^ 1 and in that special (isothermal) case we find 
P DO = F(K)e 2e , (87) 
where now 

and with v as above 
1 

— = ± . 89 

dv 2vy/2n — v — 1 00 In v 

Now it is readily found by tracing back through the various definitions that this coarse- 
grained self-similar or power law solution is a steady state. And in fact we have the relations 

Cf = (v 2 T )r-^-"' s \ C 2 2 = A- ( "- 2 " /S) , 

We can observe as a kind of verification of our procedure that in the special case a/ 5 = 3, 
the distribution function is simply / = /(«). In fact it is a member of the steady state 
solutions found for the anisotropic spherically symmetric case in ([Henriksen and Widrow 



1995| ), and used previously to model galaxies by Kulessa and Lynden-Bell ( [Kulcssa 92| ) 



although it is one of the unbound cases. It appears here in zeroth order because in this case 
as may be seen from equation fl5"TD our expansion is actually a fine graining expansion for 
a/ 5 > 2. This means that the zeroth order is the exact solution as a — > oo. For a/5 < 2 
we appear to have another set of steady solutions which develop as limits from the time 
dependent system. 

The function £(v) is not readily found analytically but it is in general an elliptic function. 
This function serves to establish a finite range in permitted values of v = £f = j 2 /r~^~ 2a ^\ 
Thus for example when a/5 = 1 we see that all large r at a given j 2 is permitted, but small 
enough r is forbidden. This is as expected. Similarly if 3/2 > a/5 > 1 so that 7 and b are 
positive but k may be negative, or if a/ 5 < 1 so that 7 is negative but b and k are positive; 
then small r is also excluded for a given j 2 . 

An obvious application of this coarse graining procedure is to the virialized core that 
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develops from a recollapsing dark matter halo. We do not now impose that the high order 
terms in the coarse graining expansion vanish. In the self-similar infall model of dark matter 
halo formation (e.g. ( [Henriksen and Widrow 1999| )) the similarity 'class' a/ 5 is given by 

a " = WtTy (90) 

where — e is the index of the density power law in the primordial density fluctuation. From 
equation (|8"2"|) we find 



dlnd _2a a (I n )X-^/a 

dh^X ~ ~T + Jloo-ilnX-WQ/aY ( ' 

Thus for example with e = 2oi a/5 = 1 and allowing the first order term Ix\X~ ( a ' s > /a 
to be as large as (1/2) J OQ , we see that the normal index of self-similar infall has flattened 
to —2 + 1 = —1. Thus the flattening at or inside this radius (where higher order terms in 
the coarse graining expansion must be taken into account)is in agreement with the NFW 
profile. However such a value for e is only reasonable on the largest halo scales. On the scale 
of galaxies the observed power spectrum index of n « —2 can be interpreted ( QHenriksen 



and Widrow 1999Q ) as requiring e = 1/2 and hence a/5 = 1/2. In this case the central 



flattening changes the self-similar index of —1 to flatter than —1/2 inside the same radius 
as above. 



5 NUMERICAL EXPERIMENTS 

In what follows we use a second generation shell code that is written in the scaled variables 
that were used in the previous sections (|Le Delliou 200 1|) , ( [Henriksen and Widrow 1999| ) to 



test the predictions of the section on radial orbits. The details of the code may be found 
in ( [Le Delliou 2001] ), but it is important to realize that these results are based in part 
(the DF figures) on an analytical estimate of the core particle potential energy in an effort 
to suppress the noise in a small number shell code simulation. This allows us to take the 
core shell potential energy to be proportional to —GM(r)/r since the additional term in 
the potential energy 4irG f£ p{r')r'dr' (where r c > is a central reference radius) is itself 
proportional to GMj r when the analytic expression for the density is used. Since the eventual 
steady state of the core appears to be reflected first in the density profiles, this seems to 
be a reasonable way to remove the noise that is due to shell crossing and arises during a 
discrete evaluation of the integral term . We are in fact able to reproduce known steady state 
distribution functions by using this technique. The second noise reducing ingredient is the 
addition of a central point mass, which assists the density profile to obtain its stable values. 
As a result of the rapid establshment of a steady core we are also able to reduce noise in the 
PDF by averaging over different stable epochs. 

We turn first to the question of whether in the radial self-similar infall phase the DF 



of a dark matter halo can really be taken as the 'one-half law' (|49|). To this end we have 



simulated the standard cosmological model of radial self-similar infall (e.g. ([Fillmore and] 



Goldreich 1984|) ; ( [Bertschinger 85|) ; (Henriksen and Widrow 1999|) ; and references therein) 



of collisionless matter using 10,000 equal mass spherical shells. We have found that the self- 
similar 'equilibrium' phase is greatly stabilized microscopically by the addition of a central 
point mass. This point mass is negligibly small (~ 5 x 10~ 4 ) compared to the mass of the 
halo so that it does not affect the global dynamics, but it imitates a central density cusp of 
zero mass as found analytically in the self-similar phase. This can not be realized numerically 
otherwise. It is particularly effective in the 'shallow' case where the halo mass is going to 
zero with the radius. 

We consider the two cases referred to as 'flat' or 'steep' in the relevant literature de- 
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pending on whether the initial cosmological density profile is flatter or steeper than r~ 2 . 
The steep cases all achieve an intermediate self-similar phase wherein the stable density 
profile is given by the initial conditions ; while the flat cases, although growing self-similarly, 
all establish a 'universal' r~ 2 density profile in the intermediate self-similar phase. That is, 
the r~ 2 profile serves as a one-sided 'attractor' for the self-similar infall. The cases that we 
show here however are selected for their stability rather than to show the evolution towards 
the attractor. This evolution has been demonstrated elsewhere at length (( [Henriksen and| 
Widrow 1999| ) and references therein). We are primarily interested here in the stable density 
and PDF profiles. 

In figure ([!]) we show the PDF averaged over the self-similar phase of the simulation 
for a 'slightly flat' case where the initial power law is oc r~ 19 . The expected core profile in 
this case is that of the attractor, namely r _2 .In the top panel we attempt to fit a negative 
temperature gaussian multiplied by an energy power law. This serves mainly as an indicator 
(at the maximum of the curve) of the limit to which the power law fit of the lower panel 
may be extended. This is the 'dominance limit' indicated on the figures. Such a law was 
suggested in flMerrit, Tremaine and Johnstone, 1989| ) and subsequently in flHenriksen and 
Widrow 1999| ) for radial infall with the power equal to 1/2. We see from the figure that this 
does not provide a very good fit to the simulated PDF. But if we focus only on the region 
where a simple power law may apply as in the lower panel, then we find good evidence for 
the expected behaviour (the slope of the averaged PDF over nearly one and one half decades 
is 0.54 ± .07). The eventual steep cut-off at large negative energies may arise in part due to 
the finite nature of the simulation wherein the most negative energies are cut off at large 
radii, but the preceding curvature may be evidence of a central (assuming that the most 
negative energies are now near the geometrical centre) gaussian PDF. We would expect such 
a cusp based on the arguments of section (Q), to the extent that the averaged behaviour 
approximates a steady state. 

In figure (|2]) the density profile established in the dark matter core near the end of the 
self-similar phase is shown. We see principally that the theoretical attractor profile of r~ 2 
is quite convincingly established over an intermediate range of scales that extends to the 
interior of the smoothing scale. We take this to be confirmation of the idea that the main 
relaxation of the core occurs near the boundary turning point since in this case more than 
the initial conditions are required to achieve the attractor profile as found. 

The profile steepens at the outside due partly to the absence of particles ( the 'Keplerian 
effect' see e.g. ( |Henriksen and Widrow 1999| )) and partly perhaps due to the fine grained 
entropy effect (the outer parts of the core retain the pure phase mixing ) suggested at the 
end of the section on radial orbits. The flattening at small radii here is largely due to the 
dominance of the point mass potential. For this purely radial model with a central point 
mass the NFW profile provides a poor fit as shown. 

In figures @ and @ the same information as above for a 'slightly steep' case is presented 
(the expected slope of 2a/ 5 = 2.03 by equation (|90|) . 

The PDF is again shown fitted by a negative temperature exponential times a power and 
by a pure power law. The power law fit yields 0.50 ± .03 when it is taken from the outer core 
regions to the maximum of the exponential curve, some one and one-half decades. However 
it is a poor detailed fit to the data here and must be considered a weaker result than is the 
shallow case. The steep cut off at high negative energies probably comes in this case from 
the finite numerical resolution of the central regions (where now the most negative energy 
particles initially, originate), but once again there is a hint of a gaussian just before this 
cut-off. 

Figure (|]) shows that the steep simulation has reproduced the expected density profile 
in the intermediate scales rather accurately, and that the profile continues well inside the 
smoothing length as anticipated. There is also the outer Keplerian steepening and an inner 



© 2001 RAS, MNRAS 000, 



Coarse Grainina the Distribution Function of Cold Dark Matter 17 



PDF with central point mass 
in shallow case £ = 1 .9 



F(E)=893.603*(-E) A 0.835887*exp(E/0.108147) 




-2 -1.5 

log(-E) 
PDF & power law 




log(-E) 

Figure 1. In both panels the data give the PDF of the simulation for a shallow initial density profile (e = 1.9) when averaged 
over the self-similar phase. The simulation was run with a small central point mass ( 5 x 1CT 4 the total mass of the halo) that 
substantially reduces the noise in the core without affecting the global behaviour. The upper panel shows an attempt to fit a 
power law times an exponential, while the lower panel fits a power law up to the point where the curvature may not be ignored. 



Shallow case £ = 1.9 

Density Profile near the end of SS phase 

io"p :.. \ 




io" 5 10 4 io~ 3 10 2 10 



Radius X 



Figure 2. Density Profile established near the end of the Self-Similar phase for an initially slightly flat perturbation. The 
smoothing scale is shown by a vertical line, but the straight line fit with slope —2 extends to the interior of this scale. The 
particles in the outer mixing zone of the halo have been neglected. The NFW profile provides a poor fit in this purely radial 
case. 
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PDF with central point mass 
in steep case £=2.1 
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PDF & power law 




Figure 3. In both panels the data give the PDF of the simulation for a steep initial density profile (e = 2.1) when averaged 
over the self-similar phase. The simulation was run with a small central point mass ( 5 x 1CT 4 the total mass of the halo) that 
substantially reduces the noise in the core without affecting the global behaviour. The upper panel shows an attempt to fit a 
power law times an exponential, while the lower panel fits a power law up to the point where the curvature may not be ignored. 



flattening due to the potential of the central point mass. Once again no globally good fit 
with the NFW profile is possible. 

We believe that these numerical experiments offer some confirmation of the ideas ex- 
pressed earlier in this paper and in (|Henriksen and Widrow 1999| ). In particular the radial 
PDF is more strongly suggested to be (— E) 1 ^ 2 in the relaxed self-similar region. This cuts 
off at large negative energies due to a combination of initial conditions and finite mass res- 
olution, but before this occurs, there is slight evidence for a Gaussian form in the central 
regions as suggested by our coarse graining in section fl2|). 

The density profiles agree with previous work in the self-similar infall phase, but we 
have reduced the noise in the simulation by the device of including a central point mass of 
negligible size. The fact that the slopes continue substantially inside the smoothing length 
suggests that the key behaviour determining the profile occurs in the mixing region near 
the boundary of the core (and v — 0) combined with the initial conditions. This also agrees 
with and supports our coarse graining argument of section It is also clear from these 
figures that the NFW profile can not provide a global fit to these radial infall models. This 
then requires an exploration of the effects of angular momentum, which we reserve for the 
next paper in this series. 
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Figure 4. Density Profile established near the end of the Self-Similar phase for an initially slightly steep perturbation. The 
smoothing scale is shown by a vertical line, but the straight line fit with slope —2.03 extends to the interior of this scale. The 
particles in the outer mixing zone of the halo have been neglected. The NFW profile can provide a good fit to the outer part 
of the halo but it is far too fiat in the inner regions for these purely radial orbits. 



6 CONCLUSIONS AND SUMMARY 

Our main interest in this paper has been to test the usefulness of a proposed coarse-graining 
(actually a progressively finer grained series) technique that uses non-canonical coordinate 
transformations on the phase space of a dynamical system. This was used in two cases in 
sections (Q) and (|3|) respectively. The first section concerned itself with the coarse graining 
of an equilibrium, spherically symmetric system by a stretching transformation on phase 
space. This allowed us to conclude that near the centre of such a system regularity of the 
coarse graining expansion imposed some conditions on the coarse-grained f{v 2 ). The one 
which maintained its form to all orders and which gave the same form at all sufficiently 
small r, was a Gaussian. We also checked that the entropy decreased in the higher orders 
(finer grained) of the expansion without imposing further constraints. The coarse graining 
expansion led naturally to the small r expansion for the density profile that is familiar from 
the theory of polytropes. 

In the second section we treated the radial self-similar infall model that is of signifi- 
cance for cosmological dark matter halos. The coordinate transformation used for the coarse 
graining was also the transformation that renders the system stationary. We found that the 
coarsest grained density profile was already that dictated by the self-similarity class implicit 
in the initial conditions. Some evidence for this behaviour was found in the numerical simu- 
lations of section (H) in that the density profile extended well inside the smoothing length of 
the simulation. We concluded that the actual relaxation of the system took place near the 
'boundary' of the core. 

The requirement that the self-similarity be stable in these coordinates to all orders of the 
coarse-graining expansion (ie essentially that the coarse grained and fine grained DF's be 
equal) yielded the self-similar DF as / oc 5(j 2 )(\E\) l l 2 . This had been previously conjectured 
in ( |Hcnriksen and Widrow 1999| ), but the present argument is more direct. Moreover from 
our simulations the evidence presented in support of this result appears stronger than before, 
although we incorporate a partly analytical estimate of the shell energy here rather than 
use the much noisier numerical calculation. Had we not done this our evidence would have 
been weaker than in QHenriksen and Widrow 1999Q . We also found that the presence of a 



point mass at the centre of the system greatly reduced the noise throughout the core in our 
simulations. 

Eventually the DF is cut off in the central regions of the system (large negative energies) 
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due to finite mass resolution and/or initial conditions. There is weak evidence for a Gaussian 
in velocity (an exponential in energy) before this cut-off, as might be expected from the 
work of secti on fl2"|). This is because the cen tral regions of the core resemble a steady state as 
discussed in flHenriksen and Widrow 1999|) , due probably to the relaxation time being short 
compared to the 'accretion time' (time to significantly change the mass) in this region. 

By expressing the entropy in the scaled variables and requiring the coarse-grained entropy 
to be positive, we saw that there was an inner limit to the self-similarity for the cosmologically 
'steep' case. Similarly there was an outer limit to the extent of the cosmologically 'flat' case. 
The limiting case which has an inverse square density profile has an exactly constant entropy 
in self-similar variables that is always positive. It diverges for an infinite system but it is 
likely to be the most stable intermediate profile for finite systems. 

By requiring the next finer grained entropy contribution to the total entropy to be nega- 
tive (corresponding to the increased information at this order) we inferred that the initially 
steep self-similar infall should flatten, primarily at the centre of the system, while the initially 
flat self-similar infall should steepen, primarily near the exterior. In this way an intermediate 
attractor behaviour exists in the self-similar infall model. However the behaviour at small 
radii is sensitive to the presence of point masses while the behaviour at large radii (exterior 
to most of the mass) tends to be Keplerian ( (|Henriksen and Widrow 1999|) ). Thus the global 
profile is more subtle in the manner of the NFW fit. However these purely radial systems 
are NOT well fitted by the NFW profile. 

We also studied theoretically (but not numerically) the spherically symmetric infall model 
in the presence of non-zero angular momentum. We showed that the coarse graining expan- 
sion can be carried out in just the same way as for radial orbits. The expected density 
profiles are not changed and flattening is expected in the central regions. The flattening is 
of a progressive nature and becomes even flatter than r _1 of the NFW profile as r — > 0. This 
is compatible with the expected Gaussian DF in the central regions of the system, although 
this can only be established definitely by examining higher orders in the expansion. There is 
an interesting dependence of the expected outer power law and flattening on scale through 
the index of the initial density perturbation (related to the power spectrum index n). This 
is such as to predict flatter cusps for galaxies than the NFW profile, which in turn is most 
appropriate on cluster scales or above (n = 1). 

We were able also to constrain the form of the DF in the relaxed state when the coarse 
grained and fine grained functions are the same as in (83) . These relaxed DF's are steady 
states as is the radial equivalent. They appear to be consistent with a possible Gaussian DF 
modulated by a function that allows for inner and outer turning points at a given j 2 . This is 
presumably due to the somewhat artificial constraint of spherical symmetry, which restricts 
the relaxation in angular momentum. 

We have found that the NFW profile can be produced in such a spherical system with 
an appropriate distribution of angular momentum (and no central black holes). However 
the angular momentum has to be correlated with initial radius of a particle in a power law 
fashion that is consistent with the self-similarity. Once again this may work only in the 
absence of strong relaxation in angular momentum. We intend to explore this question by 
coarse graining axially symmetric systems in a subsequent paper. 
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